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Q I Abstract 

I Energies of four-quark systems calculated by the static quenched 

^ ■ SU(2) lattice Monte Carlo method are analyzed in 2 x 2 bases for 

^ . square, rectangle, tilted rectangle, linear and quadrilateral geometry 

^ I configurations and in 3 x 3 bases for a non-planar geometry config- 

■ - - ' uration. For small interquark distances, a lattice effect is taken into 

account by considering perimeter dependent terms which are charac- 
terized by the cubic symmetry. It is then found that a parameter / 
- that can be identified as a gluon field overlap factor - is rather well 
described by the form exp{—[bsEA + \/57F'P]), where A and V are 
the area and perimeter mainly defined by the positions of the four 
quarks, bs is the string constant in the 2-quark potentials and E, F 
are constants. 
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1 Introduction 



QCD is thought to be the basic theory describing the interaction between 
quarks and gluons. Unfortunately, until now, because of numerical reasons, 
this can only be checked in systems containing 2 or 3 quarks. Therefore to 
make the theory more useful, methods are needed that are much less numeri- 
cal intensive than lattice Monte Carlo. These methods should agree with the 
Monte Carlo approach for the few quark systems - the assumption being that 
their natural extension to multiquark systems would also give correct results. 
Since many of the multiquark systems of interest are derived from hadron- 
hadron scattering, in which different quark clusters scatter and rearrange, it 
is the (fcf system that must be first understood. Therefore, in several at- 
tempts to describe meson-meson scatterings from QCD, the resonating group 
method in which a meson is treated as a qq cluster is a promising approach. 
In 1=2 TTTT and 1=1 KK scattering the amplitude is dominated by the direct 
gluon and quark exchange, and the quark-pair annihilation is expected to 
play a minor role. Therefore, as a first step, the detailed analysis of quark 
exchange processes from first principles or from lattice QCD is desirable. In a 
series of papers on the four-quark energies in SU(2) lattice Monte Carlo|jl|, 0], 
we considered configurations of two static qq clusters, which we called basis 
A and the quark interchanged two qq state basis B. We then analyzed the 
energy spectra in a framework similar to the resonating group method using 
the corresponding 2x2 bases. It was observed that the transition potential 
vab connecting the bases needed to be multiplied by a factor / - that can be 
interpreted as a gluon field overlap factor - in order that the solution of the 
matrix equation in this model coincided with the two lowest eigen energies 
of the Monte Carlo calculation. 

Recently the configurations of the four-quark system were extended to in- 
clude large square(LS), rectangular (R), tilted rectangle (TR), linear(L), quadri- 
lateral(Q) and nonplanar(NP) geometries^. An analysis similar to the 
showed that, except for the NP case, the gluon overlap factor /o derived from 
the ground state energy and the /i derived from the first excited state are 
close to each other, but not identical. In the NP case, the agreement was 
worse. 

For simplicity of the model, we expect theoretically /o = /i and we tenta- 
tively attribute the small difference in /o and /i as numerical uncertainties of 
the Monte Carlo calculation. In the NP case, however, the analysis in terms 
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of 2 X 2 bases could, itself, be an origin of the difficulty. 

In the SU(2) lattice Monte Carlo approach, the transition between differ- 
ent configurations occurs through multiplications of plaquettes, which char- 
acterize the gluon interaction, with the bases wave function. In the NP 
case, some plaquettes are 3-dimensional while in the other five cases they 
are purely 2-dimensional. A calculation of the glucball spectrum in lattice 
Monte Carlo shows that the 3-dimcnsional bent plaquette contains an irre- 
ducible representation of the cubic group which violates the space parity. In 
our case gluon configurations are pinned down by the position of the quarks 
and anti-quarks and the results of the glueball calculation cannot be applied 
directly, but evenso wc expect that the symmetry of the plaquettes, which 
makes the transition from basis A to B could play an important role. 

The factor / comes from the flux-tube rearrangement. In the case of TR, 
Q and NP, the flux tubes of basis B are not straight lines on the underlying 
lattice but they are a superposition of the shortest paths along the links. The 
cubic symmetry of the plaquettes that transform the basis A to the basis B of 
TR and Q have different irreducible representations of the cubic group, but 
in the NP case, there are two plaquettes that belong to the same irreducible 
representation. The excited states of the NP case could be the hybrids i.e. 
qq with gluonic excitations. We, therefore, consider two types of hybrids 
and B~ and analyze the NP data in a 3 x 3 basis using A and . 

In order to extract parameters for use in a resonating group calculation, 
the continuum limit, in which the cubic symmetry is replaced by the rotation 
symmetry, is to be considered. Thus, parameters which are speciflc to the 
cubic symmetry are not desirable. We, therefore, parametrize the factor / 
by a superposition of an area dependent term and a perimeter dependent 
term. When the area is not a plane surface, the area that we choose is not 
the sum of the area of the plaquettes but an approximate minimal surface 
that is bounded by the links. The perimeter is kept to be the length of 
links on the plaquette, since we want to respect the cubic symmetry which 
is essential for the parametrization. Wc check whether the parameters for 
different configurations like LS and R have a common parameter when the 
two configurations possess the same cubic symmetry. The success or failure 
of this approach is given by the degree to which the various configurations 
are described by the same parameters. Since the gluon overlap factor is not 
a physical quantity, it could depend on the lattice size. We expect that the 
perimeter dependence decreases as the lattice size becomes larger and in the 
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continuum limit the area term dominates. 

In sect. 2 we present the model for fitting the data of the lattice Monte 
Carlo, and the numerical results are shown in sect. 3. Finally in sect. 4 some 
conclusions are drawn. 

2 Model 

In fitting the lattice Monte Carlo we could choose three model state bases 
0- 

• A: Two links of equal length d that connect qq on 1-3 and on 2-4. 

• B: Two links that connect qq on 1-4 and on 2-3, such that the number 
of bending points is a minimum. 

• C: Two links that connect between qq and qq on 1-2 and 3-4, such 
that the number of bending points is the minimum. This basis can be 
regarded clS cl baryon anti-baryon system in SU(2). 

As explained in the Appendix of P the three states are constructed by 
contracting the operator 

Tijki = Piriz, Xi)qr{xi)Pjs{z, X2)qsix2)Pktiz, Xs)qtixs)Piuiz, Xi)quiXi)\0 > 

(1) 

with tensors 

Aijki = ^ijekh Bijki = euekj and Cijki = eikeji. (2) 

Here, in SU(2), the quark field is described by qi{x) and the antiquark field 
by Qiix) = —eijqj{x) and Pij{z,x) is the colour fiux path from z to x. Due 
to the equation eike^i = 5ij5ki — SuSkj, the base C is equal to A-B if z can be 
chosen to be a common point of the paths of B. This is the case for linear (L), 
base B2 of quadrilateral(Q) and base B]_,B2 and of non-planar(NP)(see 
figs. 4,5 and 6). If there is not a common point, Pki{z,x) for B as well as 
for A should contain a product of links that starts from z and goes to a 
quark in B and returns to z, and that starts from z and goes to a quark in 
A and returns to z, respectively. A difference in this expectation value for A 
and B would make C different from A-B, but numerical calculation indicates 
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that C is approximately expressed as a linear combination of A and B, since 
combinations A+B, A+C, B+C give very similar energy spectra|Q. In our 
colour SU(2) static model, we assign two meson states to A and B, and a 
baryon anti-baryon state to C. Therefore, base C is regarded as redundant 
and we extract the effective potential in the space A, in the space B and the 
transition potential between A and B from the lattice numerical data. 

In the actual calculation the links are fuzzed^]. Here, fuzzing means 
replacing a link Uij by a linear combination of Uij = e*^""^''" and a staple 
UikUkiUij = (.mA^a^igoA^^a^-igoA^a ^ where /i 7^ ±z/ and v ^ A. It enhances the 
overlap of the configuration and has the same effect as tadpole renormaliza- 
tion of the gluon mass. 

In the strong coupling approximation of SU(2) gauge theory, the Wilson 
matrix for a qq system is 



< iy(L,T,/) >^ {^f^'" = exp(-6,LT), (3) 

where L is the length of the link between q and g, T is the Euclidean time, 
hs is the string tension and a is the lattice constant. 

A generalization to the qqqq system, where the q and q are on the alternate 
corners of a rectangle of sides Li and L2, is[|, H 



W 



g ^-2bsL2T 

where e is a sum over transitions at time t: 

e= - g-26«L2(T-t)g-fe,LiL2g-2fesiit 

t 

In Hilbert space the transfer matrix T, which satisfies 

Z = J {dU)e-^ = TrT"" (5) 

"-0 ij (spatial) Ss"- p,t 

has the matrix element - see ref. page 103 - 

<f/'|T|f/>=exp(4^^i?eTr(f/;rif/,^.))xexp(^5:i?eTr(?7(p))), (7) 
gtao ^ ^ gia Y 
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where p means a plaquette, is the time direction lattice spacing and a is 
the space direction lattice spacing. The bare coupling constants gs and gt 
agree to the lowest order g1 = g"^ + 0{gf). 

By a group operation g, Uij transforms to [/^ = gUij and we define 
unitary operators Rij{g)\U >= \U' >. In terms of these quantities, T takes 
the form 



T = l[iT.M9)exp{^ReTrg)) x exp(^ ^ i?eTr^(p)). (8) 

The area term exp(— 6sLiL2) in eq.(4) originates from the sum of the 
spatial plaquettes as they tile the large rectangular plaquette of Li x L2. 
Therefore, the dominant part of the transition potential between the model 
space A and the model space B is written in the form VABGxp{—bsEA) where 
A is the area of the plaquette that allows transitions between the bases A 
and B, and E is to be determined by fitting the Monte Carlo lattice data. 

When the plaquette passes through 2 quarks and 2 antiquarks then the 
degeneracy of the colour configurations or the coherence of the products of 
Uij along the perimeter becomes important. The gauge transformation g 
that connects between the degenerate Z2 bases which are eigen states of the 
Cartan subalgebra t"-6°-, where for SU(2) = ia°'/2 and t°6'^ is proportional 
to ia^/29, would be expressed as a product of the form exp(t"'9'^x^) which 
is proportional to the exponential of the length of the perimeter. Thus, in 
addition to the product of the spatial tiles, which can be obtained by taking 
g = 1, we take into account the transition between different Z2 bases which 
are proportional to the length of the perimeter. 

Incorporation of the perimeter term via renormalization e.g. via the 
self energy given by the tadpole diagram is discussed also in the large N 



expansion ||10||. Therefore, we assume that the potential containing this cor- 
rection is given in a form VABGxp{—[bsEA + y^AnV]) where V is the length 
of the perimeter and Vbs is included to make An dimensionless. This form 
has been derived from the case where the four quarks are at the corners of a 
rectangle. For this the area A and the perimeter V have clear definitions. The 
further assumption is now made that, in general, A is the natural continuum 
values - independent of the underlying lattice structure. This is similar to 
describing the two-quark potential on the lattice as V{r) = —e/r + hgT + vq, 
where r is the shortest distance between the two quarks. It is known that 
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this form is good for r > 2a. However, for smaller values of r the lattized 
form 

Virij) = -ie/rij)L + b,r + Vq (9) 

is neccessary - see ref. 

For the same reasons, those four-quark configurations, where the quarks 
are within two lattice spacings, probably contain effects arising from the 
underlying lattice. We keep the lattice structure for the perimeter term V 
and take into account this effect partially. 

2.1 Irreducible representations of the cubic group 

We consider now g as an element of the cubic group or the crystallographic 
point group. It has the effect of rotating links on a plaquette and we are 
interested in group operations which connect a product of Uij to a product 
of U[j which differ only by a cyclic permutation of the links. The Wilson 
loop at a fixed time r, or the time-local operator, is characterized by L- 
tuples Om{p) = and a class formed by all the L-tuples, which 

are identical up to cyclic permutation, is denoted by [/i, • ■ ■ , fi]- The suffix 
m specifies the orientation of the plaquette and this numbering follows the 
definition of 

Here the //s stand for unit vectors in the x-, y- and z- directions: 61,62 
and 63[§]. The charge conjugation transforms C[fi, ■ ■ ■ , /l] = [— /l, ■ ■ ■ , —fi]- 
We define 01^{p) = {fL,---,fi)- usual SU(N), C-parity positive and 
negative bases are given by ± [— /l, ■ ■ ■ , — /i]- In the case of 

SU(2), there are only C-parity positive states as the bases. 

On the lattice the bases are transformed by the proper rotations A4g of 
the crystallographic point group which are divided into five classesQ; 

1. E (identity) 

2. C, = {Cs,,C^,'} (h=a,/3,7,5) 

3. CI = (h=x,y,z) 

4. C4 = {C,H,C^h'} (h=x,y,z) 

5. C2 = {C2h} (h=a,b,c,d,e,f) 

The operations on the bases are given by Aig[fi, ■ ■ ■ , /l] = [Mgfi, ■ ■ ■ , MgfL\- 
The main difference of the classification of the type of the plaquettes here as 
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compared to other lattice calculations of the glueball spectrum ||^ is that we 
fix the configuration of the links in our model space and so we restrict Mg such 
that [Mg/i, ■ ■ ■ , Mgfi] belongs to the class of [/i, ■ ■ ■ , /l] + [-/l, ■ ■ ■ , -/i]- 
Actually, the link between q and q is fuzzed and we do not know the exact 
overlap factor between the states obtained by multiplying a plaquette and 
plaquettes properly rotated which are specified by Ck- We parametrize the 
effect of the overlap of the links on the perimeter of the plaquette by a cor- 
rection to the above area law dependent term and search the parameters by 
fitting the data. 

We specify the plaquettes that allow transitions between the spacial wave 
function of A into B in the large square(LS), rectangle(R), tilted rectan- 
gle(TR), linear (L), quadrilateral(Q) and non-planar (NP) geometries as fol- 
lows. 

1. LS 

The bases and the square plaquette (sp) for the transition, which is 
specified as are shown in Fig.l. The suffix 1 corresponds 

to the assignment of ref.]^ for the number of links I = 4. 

The proper rotations that transform 0{sp)i into itself are C22, Ciz, C^J 
and that transform into 0(sp)^ are C2X, C2y, and C2b- 

2. R 

The bases and the rectangular plaquette(rp) for the transition 
which is specified as 0{rp)i are shown in Fig. 2 

The proper rotation that transforms 0{rp)i into itself is C2Z and 
that transform into 0{rp)\ are C2X and C2y 
Here it should be emphasized that, in ref.(3), the R-geometry also 
contains the four squares with sides d = 1,2,3,4. These behave 
as in case LS above. Since the LS-case has more transformation 
options than the R-case, it is expected that squares generate a 
larger overall interaction - as is observed in the lattice data. 

3. TR 

The bases and the bent plaquettes(bp) for the transition, which is 
specified as 0{bp)7, are shown in Fig. 3. The suffix 7 corresponds to 
the assignment of for the number of links / = 6. The symmetry 
property does not change for larger /. 
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The proper rotation that transforms 0{hp)-j into itself is C2C and 
that transform into 0{bp)j are C2d and C2y 

4. L 

The bases of the hnear configuration are shown in Fig.4. The tran- 
sition now occurs by the insertion of minimum plaquettes twice 
at the same position. We call it the self-energy operation. 

5. Q 

The bases and the rectangular plaquette(rp) for the transition 
which is specified as 0{rp)i are shown in Fig. 5. The transition to 
B2 occurs through the self-energy mechanism as in case L. 

The proper rotation that transforms 0{rp)i into itself is C22 and 
that transform into 0{rp)\ are C2X and C2y 

6. NP 

The bases and the plaquettes(p=sp,rp), bent plaquette(bp) and 
twisted plaquette(tp) for the transition which are specified as 
0{p)i ± 0(^)5, olbp)r ± 0{bp)i2 and 0{tp)i are shown in Fig.6. 
The suffices 7 and 12 correspond to the assignment of for 1 = 6. 
The symmetry property does not change for larger I. The proper 
rotations among bases 0{p)i and 0{p)5 and those among bases 
0{hp)j and 0{hp)i2 can be obtained by the standard method 0. 

The proper rotations that transform 0{tp)i into itself are C3/3, C^^ 
and that transform into 0{tp)\ are C2a, C2d and C2e- 
The self-energy mechanism transforms the basis A into the basis 
5"*", which is created from A by an insertion of 0{p)i + 0{p)^ and 
0{hp)'; + 0{hp)i2, but not with B~ which is created from A by an 
insertion of 0{p)i — 0{p)^ and Oihp)^ — O (bp) 12. 

2.2 The gluon overlap factor 

Among the proper rotations that keep the orientation of the plaquette, the 
operation that transforms to 0{pY, which is related to the charge conju- 
gation, would have the main overlap to the basis C. Typical examples are 
obtained by transforming the shortest qq link in the basis B2 of the NP ge- 
ometry and switching the connection of the open circle on the shortest link 
to make a qq link and a qq link. Another example can be produced in the 
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basis ^3 of NP, but this kind of coupling of the basis C to the bases A and 
B occurs only in very restricted circumstances: a diagram with three paths 
from a g or a g must be there and the three paths should be colour neutral. 
Therefore we neglect the rotation that transforms 0{p) to 0{pY. 

We define the length of the two paths in the state A and the distance 
between the two, in units of the lattice spacing a, as d and r, respectively. In 
the case of L, Q and NP, r is defined as the distance between quarks 3 and 
4. 

The strong coupling potential in the direct channel is obtained from 



<WiL,T,g') >=e 



-v{L)T 



(10) 



which yields f (r^ 



Including the lattice Coulomb interaction and the 



self-energy t>o|0 and defining Vij as (9), we parametrize va = fis + V24,vb = 
f 14 + V23. 

We fit the gluon overlap factor / for LS, R, TR, L and Q by the condition 
that the solution A of the equation 



detW = det{ 



Va 

fVAB 



fVAB 
Vb 



X 



1 

//2 



//2 
1 







(11) 



agrees with the energy eigenvalues of the lattice calculation 0. It should be 
added that the parametrization of the two body potential (9) is done for each 
configuration separately and the deviation of the lattice data from the fitted 
formula (9) is less than 0.5%. 

In the NP case, we need to analyze the lattice calculation data using the 
bases A, and . We define the overlap factors between A and -B^ by 
Since 



< 0(p), + 0ip),\C2k\0ip)i - Oip), > 

- < 0ipUC2k\0ip), > + < Oip),\C2k\0{p), > 

- < Oip)., - 0{p),\C2k\0{p), + 0{p), > 



(12) 



and the matrix W should be Hermitian, we define the transition matrix 
between B~ and is proportional to the Pauli matrix —ay multiplied by 
a factor t and choose to be complex. 

We fit and t from the condition that the solution A's of the following 
equation coincide with lowest three energy eigenvalues of the lattice calcula- 
tion. 
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V2vABf+ V2vABf- 



detW = det{ 



V2vABf+* 
V2vABf-* 




(13) 



1 



v^/+/2 V2/-/2 1 



-A 



\/2/-72 



1 



1 ) = 



(14) 



The matrix elements between A and are multiplied by ^/2 due to the 
normalization of the bases of . Algebraic forms of the three A's can be 
obtained by using Mathematica. They are functions of \f^\'^ + and 
Im{f^ f^*). Therefore and /~ are interchangeable. 

The exphcit forms of / are chosen as follows. 



Here hg is the string tension and A = (f - the area of the square. 
The parameter E fixes the overlap of the gluon configuration due 
to the strong coupling area dependence. Parameters Ai and A2 fix 
the overlap due to the links on the perimeter of the area and V — 
Ad. Ai contains contributions of and A2 contains contributions 
of C42 and C^^. 

In principle, one could include an overall normalization factor /c 
in /. In fact, if one omits the data of (r, d) = (1,1) and fits the 
overlap factor / with the additional factor /c, we obtain a much 
smaller with /c ~ 0.87. However, in this case / becomes almost 
Gaussian and the value of / at r = and r = 1 in lattice units 
becomes very close to each other. This means that when we change 
the lattice spacing a to a smaller a' the slope of / at (r, li) = (1, 1) 
and at (r, d) — (2, 2) in the new lattice spacing would be quite 
different. We compare the data of /3 = 2.4 which corresponds to 
a = 0.12/m and (3 = 2.5 which corresponds to a' = 0.082/m. The 
parameters of the quark-antiquark potential for the two f3 values 
are given in Table 1. The overlap factor / of {r,d) = (1,1) for 



1. LS 




(15) 
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(3 = 2.5 turns out to be 0.94. Thus, we fit / as an exponential and 
not as Gaussian, and fix the normahzation to be 1. This means 
that we incorporate the lattice effect or correction to the area term 
explicitly into the fitting. We then need to extrapolate the fitted 
data to the continuum limit, where the perimeter term - except 
for the self energy term - would become unimportant. 



Here A = rd and V = 2{r + d). Since R does not have the 
symmetry corresonding to C42 and C^^^, the A2 must be omitted 
here. Since the area term dominates over the perimeter term as 
r and d becomes large, the difference in the overlap factor in LS 
and R becomes small for large plaquettes. 



We define the projection of the paths 1-4 and 2-3 onto the x- 
y plane and onto the y-z plane by {x, z) - in lattice units a . 
Therefore the length of the paths in B is yx^^+z^ = r' and A' = 
dr' is the area of the surface that is bounded by the four links 
1-3, 3-2, 2-4 and 4-1. We do not take the sum of the area of 
the two plaquettes of Fig. 3 since for large d, and r it would not 
correspond to the continuum limit of interest. The spectrum for 
(d, x, z) = (5, 3, 4) = (5, 4, 3) becomes very similar to that of LS 
with r = (/ = 5, in this model. 

Although a parametrization specific to the cubic symmetry is not 
desirable, the length of the perimeter is tightly related to the cubic 
symmetry of the links: i.e. A3 contains the contribution of €20- 
Thus, we take V = 2{x -\- z -\- d) which is the length of the bent 
plaquette. 



2. R 



/ 



(16) 



3. TR 




(17) 



4. L 



— v/6jAo2r 



(18) 
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In order to simulate the self energy mechanism we take the perime- 
ter of the minimum area along the overlapping links of the basis 
B. Normalizing / = 1 for r = 0, we obtain the above formula. 

5. Q 

where the first term comes from the same mechanism as the linear 
(L) and the latter is the rectangular plaquette contribution. The 
area = rdj^'vs, the area of the triangle 1-3-4, that is half of the 
plaquette area. P = 2(r + c?) is the length of the perimeter of the 
rectangular plaquette which transforms under C2z- 

6. NP 

In the nonplanar case, combinations of bent plaquettes 0{bp)i 
(z = 1, ■ ■ ■ , 12) lead to different parity states^]. In our restricted 
plaquette bases the symmetric and antisymmetric combinations of 
the plaquettes 0{p)i ± C(p)5 and 0{hp)j ± 0{bp)i2 lead to bases 
and B~ respectively. 

We parametrize = {fre ± ifim), where 

^^-[b.EA"+Vb:A,P,]yQ ^20) 

and 

J.^ ^ ^^-[bsEAi+Vb^BiVi] _^ ^-[bsEA2+Vb:B2V2]y2 (21) 

Here Bi in fre contains the contribution of C2z ± C2y, while that 
in fim the contribution of €4^, C^^ and C2e- Vi = 2(r + d) is 
the perimeter of the plaquette 0{p)i and O^p)^, and Ai = rd 
is the corresponding area. B2 in fre contains contributions of 
C2C ± C*2b, while that in fim the contribution of C3a,C^^, and 
C2e- ^2 = 2(r + 2d) is the perimeter of the plaquette 0{hp)j and 
C(p)i2, and A2 = d{r + d) is the corresponding area. contains 
the contribution of C2c±C2f, and P3 = 2{r + 2d) is the perimeter of 
the plaquette 0{tp)i. For the corresponding area A!' we adopted 



12 



the area of a curved surface 



A" = [ du I (if |(Mri3 + (1 - M)r42) X (fr23 + (1 - t^)r4i)| 

JO JO 

= duj^ dvd^Jr^il -2u + 2u^) + d^v^ (22) 
which can be calculated analytically. 

The ansatz (22) corresponds to an approximation by the regular 
surface or an approximation of a surface by a sum of straight 
lines that connect between two flux lines in a 3-d space e.g. 23 
and 14 in fig. 6. If the four links that surrounds the area are 
not on a plane, the regular surface is not necessarily the minimal 
surface. In fact except for the plane surface it is only the helicoidal 



surface which is minimal and regular [|12|. A numerical calculation, 



using the method of |TT| shows, however, that the ansatz (22) is 



a reasonable approximation especially when r » d. Since the 
evaluation of (22) is simple, we adopt this formula for the area of 
the NP geometry. 

The appearance of t in the potential means that the configuration 
of quarks and anti-quarks are not enough to classify the energy 
spectrum and we need to specify the gluonic degrees of freedom. 
It is natural in the flux-tube picture since the links in B are bent 
and the excited state can be mixed in. The transition potential 
connects two different symmetry states which we specified by 
and B~ . 

The appearance of perimeter dependent factors and the splitting 
of B into and B~ come from gluonic degrees of freedom, since 
the lattice data are not expressed in the bases which contain only 
the quark degrees of freedom. 

3 Numerical results 

In order to fit the parameters Ai of LS, R, TR, L, and Q, we tried two 
options; 

1. Fit eigen-energies of all the {r,d) sets directly. 
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2. Calculate the overlap factor /o from the ground state eigen en- 
ergy and /i from the first excited state eigen-energy by solving a 
quadratic equation for a fixed {r,d) set|I| and fit the averages of 
/o and /i of all the (r, d) sets. 

The two options give slightly different results. The error-bars of the eigen- 
energies at small (r, d) are crucial for the relative magnitudes of the area 
dependent part and the perimeter dependent part. The error-bar for small 
{r,d) is smaller than that for large (r, d), but for {r,d) = (1,1), there are 
problems in fitting the two-body potential Vij as a sum of the three terms 
given (9)[|^. However, we assume that the lattice effect is well simulated 
by the perimeter term, and use the small error-bars given by the lattice 
simulation. The results of the first option i.e. direct eigen-energy fit give 
following parameters. 

1. LS 

We obtain E = 0.296(11) and Ai + A2 = 0.080(2). 

2. R 

We fix E=0.29Q and obtain Ai = 0.057(1) which means A2 = 
0.023(3). 

3. TR 

We again fix E=0.296 and obtain ^3 = —0.091(6). In order that 
the two mesons does not interact for large r it is necessary that 
linir^oof = 0. However, the negative value of As does not cause 
trouble since the area term dominates in physically reasonable 
configurations and / tends to for large r. 

4. L 

We obtain Ao = 0.197(3). 

5. Q 

We fix E=0.296 and Aq = 0.197 then we find A[ = 0.051(2) which 
is close to the value Ai = 0.057(1) obtained in R. This is a good 
indication that our parametrization is meaningful. When we use 
the perimeter V = d + r + \/r^~+~d? which corresponds to the 
continuum limit, we obtain for d = 2,3 and 4, A[ = 0.0595 i.e. 
closer to Ai. We insist, however, on using the V that respects the 
cubic symmetry. 
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6. NP 

In order to fit the parameters Ai and Bi of tlie NP geometry we 
first fit + = 2{fl + ffj, Imif+f-*) = 2Uehm and 

t for each set of (r, d) independently so that the three A coincide 
with the three lowest eigenenergies. In the case of > r + 1 the 
unperturbed energy of the configuration B is lower than that of 
A and we found difficulty in fitting the data. In these cases one 
should consider and B as the basis, but since the flux tubes 
of the basis B are not along the lattice link it is not possible to 
specify the bases A^. We omit these data from the analysis, except 
(r, c?) = (1,2), where the difference in the unperturbed energy is 
of the order of 0.05 and we can find parameters that fits the data 
approximately. 

We again fixed E — 0.296 and Aq — 0.197, and searched parame- 
ters to fit the I /+ 1^ + 1/- p and Im{f+f-*) of the data sets of (r, d) 
that satisfy r>d~l. We obtained Bi = 1.63(72), B2 = -0.12(1) 
and A4 = 6 ± 7. A4 consistent with zero means that the coher- 
ence of the links on the twisted plaquette is weak. A fit with /I4 
fixed to be zero gives Bi = 5.3 ± 4.7 and B2 = -0.092(13). Al- 
though the error-bar of Bi is large it is consistent with the previous 
1.63(72). The coherence on the bent plaquette parametrized by 
B2 is relatively stable. It is the main term that transforms among 
different symmetry states. We remark that B2 in the NP case 
and A3 in the TR case are both related to the bent plaquette and 
the magnitudes are similar, while Bi and Ai are both related to 
the rectangular plaquette but the magnitude of Bi is larger. In 
the latter case, not all the 4-quarks are sitting on the plaquette 
B2 and ^3. These plaquettes can be obtained by multiplying a 
plaquette to the base B4, and by multiplying a plaquette to the 
base B5, respectively, on which all the 4-quarks are sitting. Since 
a link in the model bases contains a staple-shaped link, due to the 
fuzzing, the effective area and the length of the perimeter of the 
plaquette characterized by the parameter B2 would be modified. 
We cannot, however, estimate the correction and the parameter 
Bi is regarded as a phenomenological parameter. 
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The parameter t that connects the bases and -B" should tend 
to as r and d increase, and we observed this tendency in the 
fitting. But for simphcity, we search an average value of the pa- 
rameter t that fits the three lowest eigcncncrgics of all the data 
of {r,d). Taking Bi = 1.63,52 = -0.12 and A4 = 6.19 and the 
error-bars of the cigcn-cncrgics given by the lattice calculation 
data, we fitted t. Due to relatively large error-bars in the highest 
energy state, we sacrifice the accuracy of the fit of the highest 
state. In fact in the highest state, more complicated configura- 
tions could be involved and the fitting with three bases could be 
a bad approximation. The value we obtained for t is 0.045(10). 

4 Discussion and Conclusion 

We analyzed the four-quark energies in SU(2) calculated by lattice Monte 
Carlo in a model which is inspired by the resonating group calculation of 
hadron interactions. The transition matrix element between different bases 
A and B is assumed to be proportional to the gluon overlap factor and this 
factor is then parametrized by a product of an area dependent and a perime- 
ter dependent term. The perimeter dependent term was introduced by con- 
sidering transition between degenerate Z2 symmetric bases which in the large 
lattice or in the continuum limit would vanish except the self energy term. 
In the case of large squarc(LS), rectangular (R), tilted rectangular (TR), lin- 
ear(L) and quadrilateral (Q) configurations, we found reasonable parameters 
for the perimeter dependent term which are characterized by the cubic sym- 
metry. In the non-planar (NP) case, we observed that it is necessary to take 
into acount the gluonic degrees of freedom explicitly. The base B was split 
into two hybrid configurations: and B~ . Using A and B^ as the bases, 
we parametrized the gluon overlap factor for the NP. 

We treated, the error-bars given by the lattice Monte Carlo as meaningful 
also for the small plaquettes and checked whether a consistent parametriza- 
tion for different geometries is possible. We simulated the deviation from an 
area dependence by adding a perimeter dependent contribution. 

The parameters for the perimeter dependent part are Aq, Ai, A2, A^, A4, Bi 
and B2. A[ in the Q case is consistent with that of R case, and B2 in the 
NP case is consistent with A3 in TR case. These coincidences are not trivial 
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and are encouraging for any further analysis based on the cubic symmetry 
of the flux tube geometry. A4 is consistent with and we could well omit it. 
Therefore, there are 5 parameters for the perimeter terms and one parameter 
for the area to fit the lattice Monte Carlo energy spectra. 

We expect that in the continuum limit / is completely given by the area 
term except for the self energy contribution. A comparison of the data at 
P = 2.4 and /? = 2.5 for the LS configuration suggests that the parameters of 
the perimeter term in fact decrease as the lattice constant a becomes small. 
Since the Monte Carlo data contain lattice artefacts, in order to analyze the 
data for a specific configuration we need to choose proper bases so that these 
artefacts are eliminated. Since the artefacts are represented by the perimeter 
term, its cubic symmetry property is essential for choosing the proper bases. 
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Figure 1: The basic four-quark geometries of Large square(LS). A is the 
configuration of the base A. Bi are the configurations of the base B. The 
hnked qq of the bases A are assigned 1-3 and 2-4, and those of the bases B 
are 1-4 and 2-3. 



Figure 2: Same as Fig.l but for Rectangle(R). 



Figure 3: Same as Fig.l but for Tilted rectangle(TR). 



Figure 4: Same as Fig.l but for Linear(L). 



Figure 5: Same as Fig.l but for Quadrilateral(Q). 



Figure 6: Same as Fig.l but for Non- planar (NP). 



Table 1: bs,e and Vo of the quark-antiquark potential 





bs 


e 


"^0 


2.4 
2.5 


0.07169 
0.03781 


0.245 
0.215 


0.550 
0.523 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9409006vl 



